Critical charge instability on verge of the Mott transition 
and the origin of quantum protection in high-Tc cuprates 



o 
o 

(N 



o 



> 

(N 

cn 
cn 

o 
o 



I 

o 
o 



X 
J3 



T. K. Kopec 
Institute for Low Temperature and Structure Research, 
Polish Academy of Sciences, 
POB 1410, 50-950 Wroclaw 2, Poland 
(Dated: February 6, 2008) 

The concept of topological excitations and the related ground state degeneracy are employed 
to establish an effective theory of the superconducting state evolving from the Mott insulator for 
high-Tc cuprates. The theory includes the effects of the relevant energy scales with the emphasis 
on the Coulomb interaction U governed by the electromagnetic U(l) compact group. The results 
are obtained for the layered t — t' — t±— U — J system of strongly correlated electrons relevant 
for cuprates. Casting the Coulomb interaction in terms of composite-fermions via the gauge flux 
attachment facility, we show that instanton events in the Matsubara "imaginary time" , labeled by a 
topological winding numbers, are essential configurations of the phase field dual to the charge. This 
provides a non-perturbative concept of the topological quantization and displays the significance of 
discrete topological sectors in the theory governed by the global characteristics of the phase field. 
We show that for topologically ordered states these quantum numbers play the role of an order 
parameter in a way similar to the Landau order parameter for conventionally ordered states. In 
analogy to the usual phase transitions that are characterized by a sudden change of the symmetry, 
the topological phase transitions are governed by a discontinuous change of the topological numbers 
signalled by the divergence of the zero-temperature topological susceptibility. This defines a novel 
type quantum criticality ruled by topologically conserved numbers rather than the Landau principle 
of the symmetry breaking. We show that in the limit of strong correlations topological charge is 
linked to the average electronic filling number and the topological susceptibility to the electronic 
compressibility of the system. We exploit the impact of these nontrivial U(l) instanton phase field 
configurations for the cuprate phase diagram which displays the "hidden" quantum critical point 
covered by the superconducting lobe in addition to a sharp crossover between a compressible normal 
"strange metal" state and a region characterized by a vanishing compressibility, which marks the 
Mott insulator. Finally, we argue that the existence of robust quantum numbers explains the stability 
against small perturbation of the system and attributes to the topological "quantum protectorate" 
as observed in strongly correlated systems. 



I. INTRODUCTION 



Superconductivity represents a remarkable phe- 
nomenon where quantum coherence effects appear at 
macroscopic scakfi The superconducting (SC) proper- 
ties, especially the perfect diamagnetism, are microscopic 
manifestations of the spontaneous breakdown of one of 
the fundamental symmetries of matter, namely the U(l) 
gauge symmetry.'^ For a superconductor this introduces 
a state with no definite particle number, but with a def- 
inite conjugate variable identifiable with the phase. The 
discovery of the cuprate superconductors,"^ which is be- 
lieved to emerge due to the strong electron-electron in- 
teractions, has sparked a widespread interest in physics 
which goes beyond the traditional Fermi-liquid frame- 
work usually employed for understanding the effect of 
interactions in metals>J' It has been recognized that the 
superconductivity occurs in the region of the doped Mott 
insulator near the Mott transition, so that the whole mi- 
croscopic foundations on which BCS theory was built on 
cannot be accommodated to explain cuprate supercon- 
ductivity. The conventional BCS theory identifies SC or- 
der as an instability of the Fermi sea. That assumes the 
existence of a normal Landau-Fermi-liquid^ which forms 
well above the critical temperature Tc, which eventually 



turns into a superconductor below Tc due to residual at- 
traction among low-lying quasiparticles. In cuprates, in 
turn, there is overwhelming evidence that superconduc- 
tivity does not appear as an instability of a Landau- Fermi 
liquid. The reason is that the fermion quasiparticle does 
not seem to reflect the character of the measured low en- 
ergy eigenstates. In this sense, the electron (or electron- 
like quasiparticle) may no longer be the appropriate way 
to think of elementary excitations. It is intriguing to 
conjecture that in the strongly correlated systems the 
ubiquitous competitions between the variety of possible 
ground states govern the essential physics - the forma- 
tion of a highly degenerate state seems to open the way 
to transformation into alternative stable electronic con- 
figurations. Strong correlations that suppress electron 
motion, may transform the system into a kind of a un- 
stable state which will be very sensible to charge and/or 
spin ordering. In this critical-like state the superconduc- 
tivity might emerge as a competition between different 
ground states. Indeed, in cuprates there is clear evidence 
for the existence of a special doping point "hidden" near 
optimal doping below the SC dome where superconduc- 
tivity is most robust Such a behavior suggests that 
this point could be a quantum critical point (QCP) while 
the associated critical fluctuations might be responsible 
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for the unconventional normal state properties»& How- 
ever, it is unclear whether this QCP is "truly critical" 
in the sense that it is characterized by universality and 
scaling. For example, if excitations at QCP also carry 
the electrical current, then a resistivity linear in tem- 
perature, as is observed in cuprates, obtains only if the 
dynamical exponent is unphysically negatived The re- 
semblance to a conventional QCP is also hampered by 
the lack of any clear signature of thermodynamic criti- 
cal behavior. Usually, a QCP would be the end-point of 
a critical line below which an ordered phase takes place 
and it could be made manifest below the superconduct- 
ing dome. Experiments appear to exclude any broken 
symmetry around this point although a sharp change 
in transport properties is observed.^" Unfortunately, our 
understanding of the underlying orders in cuprates is far 
from being satisfactory and identifying the nature of the 
putative QCP is an open question. A possible way out 
from this difhculty would be if the degrees of freedom that 
rule QCP's are different from the energy degrees of free- 
dom that govern the stable phases the critical point sepa- 
rates. Thus these "non-conventional" degrees of freedom 
could provide critical fluctuations beyond those of the or- 
der parameter fluctuations usually included in the stan- 
dard Ginsburg-Landau- Wilson (LOW) descriptioui ^^i^^ 
Within the LGW approach the order parameter fluctu- 
ations arc implemented by constructing models which 
mimic the low energy properties of solids - a procedure 
which relies on the separability of the high and low energy 
scales. Under this proviso the high energy variables can 
be removed out yielding an effective Hamiltonian (as ex- 
emplified e.g. by the t — J model^^) which describes the 
relevant low energy and long- wavelength physics. This 
procedure is often implemented via the projective trans- 
formation, which results in removing of high-energy de- 
grees of freedom and replacing them with kinematical 
constraints. In such approaches, the high energy scale 
associated with the charge gap is argued to be irrelevant, 
hence the focus exclusively on the spin sector to char- 
acterise the Mott insulator. However, the charge trans- 
fer nature of the cuprates plays an essential role in the 
doped systemsii so that with discarding charge degrees 
of freedom an important part of the physics may be lost. 
However, there is also mounting experimental evidence 
which put in question the validity of the various projec- 
tion schemes. For example, it has been found that above 
any temperature associated with ordering in both elec- 
tron and hole-doped cuprates a charge gap of order 2eV is 
present in the optical conductivity and a rapid reshuffling 
of spectral weight with hole doping. ^^■-'^^ Angle-resolved 
photoemission spectroscopy (ARPES) also reveals a sim- 
ilar charge ga,pM. Surprisingly, when superconductivity 
emerges, the low and high energy degrees of freedom are 
still coupled. It has been shown that changes in the op- 
tical conductivity occur at energies 3eV which exceeds 
by two orders of magnitude the maximum of the super- 
conducting gap,"^ But the high energy scale involved is 
hard to understand unless it is assumed that the spec- 



tral weight is transferred from both the lower and the 
upper Hubbard bands, thus beyond the range of appli- 
cability of the pure t — J model. The intimate connect- 
edness between the low and high energy degrees of free- 
dom in doped Mott insulators was firmly appreciated and 
termed mottness^^ Due to the intrinsic mixing of high 
and low energy degrees of freedom no low energy reduc- 
tion is possible in a conventional sense, so that the doped 
Mott insulators are inherently asymptotically slaved.^" In 
a similar spirit a detour from the strict projection pro- 
gram was recently proposed in a form of the "gossamer" 
superconductor^" recognizing the role of the expensive, 
double-occupancy charge configurations. 

While spontaneous symmetry breaking has become one 
of the main guiding principles in physics there are 
other signatures in a physical system that are associated 
with the topological effects. These are instrumental for 
a full understanding of the physics and lead to a host 
of rather unexpected and exotic phenomena, which are 
in general of a nonperturbative nature. For example, 
the fundamental character of a vector potential is ev- 
ident in the Aharonov-Bohm (AB) eSect^^ where the 
topology of the U(l) group is essential: when an elec- 
tron is transported in a magnetic field around a closed 
loop, it acquires a phase that is equal to the magnetic 
flux through the surface spanned by the electron path. 
Strongly correlated electronic system are no exception in 
this regard. In particular, the fractional quantum Hall 
(FQH) effectSiiS^ is the prominent representative. Here, 
the striking fact is that FQH systems may contain many 
different phases at zero temperature which have the same 
symmetry. Thus different states cannot be distinguished 
by symmetries and the Landau symmetry-breaking prin- 
ciple fails because also topological characteristics of the 
configurational space come into play. In the most inter- 
esting cases configurational spaces are not simply con- 
nected. There are space-time configurations of quantum 
fields which cannot be continuously deformed one into 
another. Further, an adiabatic motion along a noncon- 
tractible closed path in a configurational space leads to a 
geometric (Berry) phase^^ acquired by the wave function. 
In condensed matter systems with large numbers of mu- 
tually interacting particles, the subject of Berrys phases 
becomes a key issue. For example, a Berry phase distin- 
guishes between integer and half-integer spin chains and 
results in different ground states and excitations,^^ The 
other examples are the possible geometric phases effects 
on statistical transmutation^^ that can be achieved by 
a "fiux attachment" which now becomes a very power- 
ful theoretical method. In many cases the topological 
character of the quantum field is captured by single inte- 
ger, called the topological charge, or winding number of 
the field which classify topological excitations. These are 
found by integrating the so-called Chern-Simons terms 
which enter the Lagrangians of the theory. What is com- 
mon to all the above issues is the appearance of gauge 
fields to characterize various interactions: field configu- 
rations which differ by a gauge transformation are to be 
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regard as physically the same. 

In the present work we argue that the important prop- 
erties of cuprates are controlled by the large Mott gap 
and consider the representation of strongly correlated 
electrons as fermions with attached "flux tubes". This 
introduces a conjugate U(l) phase variable, which ac- 
quires dynamic significance from the electron-electron 
interaction. This means that an electron is not a quasi- 
particle (in the Landau sense), but has a composite na- 
ture governed by the electromagnetic gauge group. Fur- 
thermore, we recognize the non-trivial topology of the 
electromagnetic U(l) group by observing that the funda- 
mental group 7ri(?7(l)) = Z is given by a set of integers. 
Therefore the elementary excitations in a strongly cor- 
related system always carry 27r-kinks of the phase field 
characterized by the topological winding number^''^ - a 
quantized U(l) topological charge. Due to the nontrivial 
first homotopy group tti the configuration space of the 
quantum phase fields is multiply connected, so that in- 
equivalent paths in the "imaginary-time" (paths that are 
not deformable to one another) naturally emerge. Hence, 
the U(l) gauge field gives rise to a "topological interac- 
tion", which is felt by the electrons and can be sepa- 
rated from ordinary dynamical ones lumping them into 
particle "statistics". To facilitate this task we employ 
the functional integral formulation of the theory that en- 
compasses all of these topological possibilities: one has 
to perform the functional integral over fields defined on 
different topologically equivalent classes i.e, with differ- 
ent winding numbers. From a canonical (operator) point 
of view, however, the different topological sectors seem 
to give rise to completely different Hilbert spaces and 
the resultant field operators would satisfy quite compli- 
cated nonlocal commutation relations. The fact that a 
prospective theory of electronic states in strongly corre- 
lated electron systems must give up on either standard 
fermion commutation relations or standard particle con- 
servation laws has already been pointed out.^^ Further- 
more, we exploit the impact of these topological excita- 
tions for the phase diagram of cuprates (with its various 
crossovers and transition lines) and show that they can 
induce its unusual feature: a "hidden" quantum criti- 
cal point of the type that results from the topological 
ground-state degeneracy. It can be probed by the topo- 
logical susceptibility as a robust, nonperturbative prop- 
erty that is related to the physical quantity of interest, 
namely, the diverging charge compressibility. It also pro- 
vides natural description of the Mott state where the sys- 
tem is said to be incompressible when there is a gap in 
the chemical potential as a function of the electron den- 
sity. This topological underpinning establishes the source 
of quantum protection as a collective state of the quan- 
tum field, whose excitations pertain to the whole system. 
Therefore, macroscopic behaviour is mostly determined 
by topological conservation laws which does not arise just 
out of a symmetry of the theory (as "conventional" con- 
servation laws based on Noether's theorem) but it is a 
consequence of the connectedness, i.e. topology of the 



phase space, related to the topological properties of the 
associated gauge group manifold. 

The organization of this paper is as follows. In Section 
II we introduce the electronic model for cuprates which 
captures the layered structure and epitomizes the hier- 
archy of relevant energy scales, with the largest set up 
by the Coulomb interaction. In Section III we describe 
the details of the flux attachment transformation which 
results in the representation of strongly correlated elec- 
trons as fermions plus attached U(l) gauge "flux tubes" 
that leads to composite particle picture. Section IV is 
devoted to the basic concepts of the algebraic topology 
(homotopy groups) that becomes instrumental for the 
Feynman's path integral formulation of quantum statis- 
tics. This is followed by Sections V and VI where the 
fermionic part of the theory is elaborated in terms of 
the momentum dependent "d-wave" spin-gap and mi- 
croscopic phase stiffnesses. Here, the most important re- 
sults are summarized in the effective bosonic model writ- 
ten with the help of the collective phase variables. This 
enables us to study the superconductivity as the conden- 
sation of the "flux-tubes" from the electron composite 
which is presented in Section VII. We elucidate there the 
role of doping for the superconducting order to occur and 
the key role played by the topological degeneracy. In the 
subsequent Section VIII the topological susceptibility is 
used to probe the change of the topological order. Here, 
we show that the topological susceptibility is related to 
the charge compressibility that diverges at the degener- 
acy point at zero-temperature and defines a novel type 
of topological quantum criticality, beyond the paradigm 
of the symmetry breaking. In Section IX we present cal- 
culated phase diagrams for cuprates displaying, beyond 
the conventional ordered states, regions that are related 
to the change of the topological order. Section X is ded- 
icated to the discussion of the robustness of the ground 
states in cuprates and its source in the topological con- 
servation laws. Finally, in Section XI we conclude with 
a summary of the results, while the Appendixces collect 
material that is related to the technical part of the work. 



II. ELECTRONIC MODEL FOR CUPRATES 

The Hubbard model is viewed as the generic model 
for interacting electrons in the narrow-band and strongly 
correlated systems^ that captures the physics of Mott 
transitioni^ The existence of the Mott insulator in the 
cuprates' parent compounds implies that a viable theory 
of high-temperature superconductivity must explicitly 
incorporate the Mott-Hubbard gap for charge transfer ii4 
While it is believed that the basic pairing mechanism 
in cuprates arises from the antifferomagnetic (AF) ex- 
change correlations^ it is apparent the charge fluctua- 
tions also play an essential role in doped systems. Hence 
the Coulomb charge fluctuations associated with double 
occupancy of a site are controlled by the parameter U in 
the Hubbard model, which also determines the strength 
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of the AF exchange couphng J. Therefore energy scale 
of the charge fluctuation is characterized by the Mott 
gap, which is by far larger than the energy scale of mag- 
netic fluctuations. Although the t — J model is usually 
viewed as the U ^ oo limit of the single band Hubbard 
model, the one-particle spectra of the two models differ 
considerably^ As already mentioned, this limitation of 
the t — J model comes from having projected out the 
doubly occupied states originally contained in the Hub- 
bard model. To explore the more flexible arrangements 
between J and U than encoded either in the Hubbard or 
t — J models we employ in the present paper a general- 
ized t — t' — U—J model for the Cu02 plane in high— Tc 
superconducting. 

In this way we retain basic features of both models: 
the charge fluctuations present in the Hubbard model 
(but removed from the constrained t — J model) and the 
robust superconducting correlations described by the ex- 
change interaction J. Furthermore, we incorporate be- 
sides the direct nearest neighbor hopping t also the next 
nearest neighbor t' hopping parameter, which importance 
in cuprates has been emphasizedi^i A good deal of the ex- 
isting literature on the cuprates invokes model Hamilto- 
nians based only on the properties of a single CuO layer. 
Obviously, the interlayer structure cannot be ignored!^ 
the measured critical temperature Tc is strongly depen- 
dent on the interlayer structure, Therefore, three di- 
mensional (3D) coupling of planes must play an impor- 
tant role in the onset of superconductivityjiMi which 
have to be incorporated by means of the interlayer cou- 
plings and the associated c-axis dispersion effects in mod- 
ehng of the cuprates^^ 

Summing up, we consider an effective one-band elec- 
tronic Hamiltonian on a tetragonal lattice that empha- 
sizes strong anisotropy and the presence of a layered 
Cu02 stacking sequence in cuprates: H = Ht+Hj+Hu, 
where 



(rr') 



i'4(r)c„,(r') - ^tx(rr')4(r)c„,+i(r') 

({rr')) rr' 



e (rr') 



SKr) • S,(r') - -m{r)ne{r') 



(1) 



Here (r, r') and ((r,r')) identify summations over the 
nearest-neighbor and next-nearest-neighbor sites la- 
belled by 1 < r < iV within the CuO plane, respec- 
tively, while I < £ < N± labels copper-oxide layers. 
Subsequently, t and t' are the bare hopping integrals, 
while t± stands for the interlayer coupling. The oper- 
ator c]j^(r)(cc(^(r)) creates (annihilates) an electron with 
spin a at the lattice site (r,£) and n^(r) = n^£{r)+nii{r) 



stand for number operators, where nQ,f(r) = cJj^(r)cQ^(r) 
and fi is the chemical potential. Furthermore, 



Q/3 



(2) 



denotes the vector spin operator (a ~ x, y, z) with ctJ^^ 
being the Pauli matrices. Finally, U is the on-site repul- 
sion Coulomb energy and J the AF exchange. Owing the 
lattice arrangement the full electronic dispersion is given 
by 

e(k,fc,) = e||(k)+ex(k,A:,), (3) 

where the in-plane contribution reads 

e||(k) = — 2t [cos(afca;) -f cos(afcj,)] 

-I- 4:t' coB{akx) cos{aky) (4) 

with t' > 0. Furthermore, the c— axis dispersion is given 

by 

e±(k, A:^) = 2i^(k) cos (c/c^), 

tj_(k) = t±_ [cos(afc2;) — cos(afcj,)]^ (5) 

as predicted on the basis of band calculationsji^ 

III. ELECTRON AS A COMPOSITE OBJECT 

We now provide the representation of interacting elec- 
trons as fermions plus attached "flux tubes" .^^ This leads 
to a picture of composite particles which are void of 
the mutual interactions among fermions: the electron- 
electron Coulomb interaction will be transformed into 
the action of U(l) gauge (phase) fields governed by the 
effective kinetic term of "free" quantum rotors. 

A. Fermionic action 

The partition function for the system governed by the 
Hamiltonian in Eq. |^ can be represented as a path inte- 
gral using fermionic coherent-states. Introducing Grass- 
mann fields Ca^(rT), CQ£(rT) that depend on the "imagi- 
nary time" < T < /9 = I/fcsT, with T being the tem- 
perature, we write the path integral for the statistical 
sum Z as 



Z 



[VcDc 



-S\c,c\ 



with the fermionic action 



f 

S[c, c] = dr 
Jo 



^ Cae{rT)drCae{rT) + n{T) 



(6) 



(7) 



The Hubbard term in Eq.JQ) we write in a SU(2) invari- 
ant way as 

Huir) = [/^{(l/4)n/(rT) - [f2,(rr) • Sdrrf} (8) 
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singling out the charge-U(l) and spin-SU(2)/U(l) sec- 
tors, where the unit vector flilrr) labels varying in 
space-time spin quantization axis.**"* The spin-rotation 
invariance one can made explicit by performing angu- 
lar integration over O(rT) at each site and time. In 
the following we fix our attention on the U(l) invariant 
charge sector leaving aside possible magnetic orderings 
such as antiferromagnetism. Although sometimes con- 
current magnetic transitions occur with the Mott transi- 
tion, the mechanism of the Mott transition is primarily 
independent of the symmetry breaking of spins. Thus we 
stress our primary interest in cuprates due to their su- 
perconducting properties and the fact that the supercon- 
ductivity (resulting from condensation of charge) should 
not be viewed as inextricably connected with the quan- 
tum antiferromagnetism. A clear support for this point 
comes from the observation of the spectral weight trans- 
fer through the superconducting transition in cuprates, 
which cannot be explained by invoking the antiferromag- 
netic order: the spectral weight transfer persists well 
above the Neel temperature and at the doping level where 
antiferromagnetism is absent. 



B. Gauge flux attachment transformation 

To proceed, we employ the Hubbard-Stratonovich 
transformation to decouple the Coulomb term in Eq.® 
with the help of the fluctuating imaginary electrochem- 
ical potential iVe{rT) conjugate to the charge number 
ni{rT). Furthermore, we write the field Veirr) as a sum 
of a static Voi (r) and periodic function^ 



Viirr) = Vi{rT + l3) 
VivT) = Vo{r) + VirT), 

where using Fourier series 

oo 

f(rr) = (l//3)^[^(rc.„)e^'^"^ + 



(9) 



(10) 



with Un = 27rn//3 (n = 0, ±1,±2) being the (Bose) Mat- 
subara frequencies. Now, we introduce the phase (or 
"flux") field (/)£(rr) via the Faraday-type relation 



?(rT) 



?(rT) 



dr 



(11) 



to remove the imaginary term 



. f'^ ~ 

I / dT4>i,{vT)ni{rT) = i drV^ (rT)nf (rr) (12) 











for all the Fourier modes of the V£{rT) field, except for 
the zero frequency by performing the local gauge trans- 
formation to the new fermionic variables faii^T): 



Caeirr) = exp 



i / dr' Viirr') 



Thus as a result of Coulomb correlations the electron 
acquires a phase shift similar to that in the electric (i.e. 
scalar) AB effect The expression in Ea. (|13|l means that 
an electron has a composite nature made of the fermionic 
part fai{YT) with the attached "flux" (or AB phase) 
exp[i(/)£(rT)]. Here, the quantity Co(rT) defined by 



Cffl(rT) = cpt{YT) = F,(rr), 



(14) 



(13) 



is the one dimensional (temporal) component Chern- 
Simons teTwM> that makes the minimal coupling with the 
fermion density field. Since the abelian Chern-Simons 
term is just as a total (time) derivative, the integral 
of it becomes simply converted by into a "surface" inte- 
gral, sensitive only to the global properties of the U(l) 
gauge field along a "imaginary time" path that starts at 
imaginary time r = and ends at r = /3. Thus the 
paths can be divided into topologically distinct classes, 
characterized by a winding number defined as the net 
number of times the world line wraps around the system 
in the "imaginary time" direction. As we shall see in the 
next Section, our considerations related to the "imagi- 
nary time" boundary conditions can be formalized us- 
ing homotopically non-trivial gauge transformations for 
which the strength of the phase shift must be quantized, 
so that the gauge change of the Chern-Simons term will 
be an integral multiple of 2tt. 



IV. QUANTUM STATES ON MULTIPLY 
CONNECTED SPACES 

A. Homotopy theory 

The algebraic topology and precisely the concept of 
homotopy groups^ provides the necessary background 
by making reference to the topological structure of the 
group manifold, let us say M . The n— th homotopy group 
7r„(M) is a group of equivalent classes of loops that can 
be smoothly deformed into each other without leaving M, 
where n refers to the dimensionality of the loops in ques- 
tion. The already mentioned AB effect is of topological 
nature, since mappings from the electron's configuration 
space to the gauge group constitute the non-trivial ho- 
motopy group 7ri(J7(l)) = Z, where the elements of Z are 
integers and represent winding numbers, i.e, U(l) topo- 
logical charges. We are precisely facing a similar situation 
in a strongly correlated system since the electromagnetic 
U(l) gauge group governs the charge sector. This be- 
comes apparent in an electrodynamic dual description of 
the charge in terms of the phase field, see Eq. Hll|) , having 
its roots in the quantum mechanical complementarity of 
phase and number. In the "imaginary time" evolution of 
the phase, two field configurations lie in the same con- 
nected component of configurational space if they can 
be continuously deformed into each other. There is a 
natural equivalence relation between these paths called 
homotopy: two paths are equivalent (i.e. belong to the 
same class) if they can be "smoothly deformed" into each 
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e(T)+2imiT/p 




FIG. 1: Schematic representation of tlie mapping from the 
real line R (the covering group of (7(1)) to the circle Si ~ U{1) 
which is locally invertible provided the topological sector with 
the winding integer number m is chosen. The map is defined 
by a continuous function f{6) on [0,27r], where /(O) = and 
continuity of the map requires that /(27r) = 27rm, so that the 
homotopy classes of 'ki{Si) are labelled by integers. 

other. The classes are labeled by the winding numbers 
and are endowed with a group structure by appropriately 
defining the composition of two mappings 4i The rule of 
thumb is that if the homotopy group is trivial, then there 
cannot be any topological field configurations in the un- 
derlying theory. 

B. Homotopy classes and the path integral 

If we work in the Feynman's path integral formula- 
tion of the quantum statisticsiS, then the statistical sum 
Z takes a form, in which homotopically distinct paths 
have to be summed according to various possibilities for 
inequivalent quantizations (superselection sectors). Spe- 
cializing to the case of the U(l) group one obtains 

^ p{m)Z{m) (15) 

me7ri((7(l)) 

Here, m £ Z labels equivalence classes of homotopically 
connected paths and p{m) marks the "statistical " weight 
which is related to a homotopy class. Furthermore, the 
partial sum Z{m) within the m-th topological sector is 



given by the usual path integral 

2(m) [ [^'^1™ [^/^/] e"""^'^"'^' (16) 

with the integration restricted to the m-th homotopy 
class. Furthermore, the weights factors p{m) form uni- 
tary irreducible representations of the homotopy group 
7Ti{U{l)), SO that the conditions for the weights are 

|p(m)| = 1 

p{mi)p{m2) ^ p{mi-km2). (17) 

In Ea. (|17|l mi and m2 label the homotopy classes of two 
the paths with common end points, while mi ★ m2 la- 
bel the homotopy class of the path obtained by joining 
the two. In particular, the weight factor, which furnishes 
the representations of Z takes the form p{m) = e'^™, 
where S [0, 27r) is the "statistical angle" parameter42. 
From the canonical point of view, if the configuration 
space of the system is not simply connected, as for the 
U(l) group, then the quantization prescription becomes 
ambiguous since paths belonging to different homotopy 
classes can get the extra relative Berry phase factor ac- 
quired by the wave function and the 9 factor represents 
exactly this quantization ambiguity. Moreover, the 9 
term cannot be traced in a perturbation theory because 
it has no affect upon the equations of motion. Therefore, 
in performing the integration in Eq. I|16|l one should take 
phase configurations that satisfy the boundary condition 
0f(r/3) - 0£(rO) = 27rm£(r) and m^(r) = 0,±1,±2, .... 
For the sake of convenience it is desirable to "untwist" 
the boundary condition by setting (see, Fig^ 

2nT 

(j>i{rT) — > <j>i{rT) EE ipeirr) + — — m£(r). (18) 

P 

Summing over all the phases (j)i{vT) amounts to integrat- 
ing over the /3-periodic field (pt,{YT) ((^^(rO) = (/?£(r/3)) 
and to sum over the set of integer winding numbers 
{mf(r)}. The integral over the static (zero frequency) 
part of the fluctuating electrochemical potential (Vo(r)) 
we calculate by the saddle point method to give 




Explicitly, following the prescription given in Eas.H15|) 
and \i<6\ we obtain for the statistical sum 



^ = E / n^^"^W / \[d^i{vT) [Vjvf]e-'^^^'"^'f'f\ (20) 
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with the action involving the topological Chern-Simons term and the statistical angle parameter 



2/1 ^ 1 
It 2^7 



Here, TC [0, /, /] is the effective Hamiltonian that is void of the Coulomb interaction 

e ^ I ra ra 



((r,r')> 



(r,r'> a (rr') 



(21) 



(22) 



where fi = ^ — nfU/2 is the shifted chemical potential, 
while 



E(/"M/"Hrr)) 



(23) 



is the occupation number for the fermionic part of the 
electron composite. In Ea. H22|l . while writing the term 
that governs AF interaction we made use of the following 
representation 

g,(rT,r'r) = -L[/^,(rr)/7,(r'r) - ^(rr)^(r'r)] (24) 

which is just the singlet -pair (valence bond) operator^'^ 
emerging from the decomposition 



S,(rr) • S,(r'r) - in,(rr)n,(r'T) 



-j5]5]e,(rr,r'T)6,(rr,r'r). 

e {rr'> 



(25) 



It is obvious that a quasiparticle description (of any kind) 
makes sense only when the constituent objects are weakly 
interacting. The chief merit of the transformation in 
Ea. (|13|l is that we have managed to cast the strongly 
correlated problem into a system of weakly interacting f- 
fermions with residual interaction given by J, submerged 
in the bath of strongly fluctuating U(l) gauge potentials 
(on the high energy scale set by U) minimally coupled to 
/-fermions via "dynamical Peierls" phase factors. It is 
clear that the action of these phase factors "frustrates" 
the motion of the fermionic subsystem. However, as we 
demonstrate in the following, it is only when charge fluc- 
tuations become phase coherent the frustration of the ki- 
netic energy is released. From Ea. (|21|l we can read off 
the explicit expression for the statistical weight 



p{m) 



exp 



« • — • 2TTme(r) 



(26) 



so that the statistical angle reads 

e _2fi 
2^^ ~U' 



(27) 



The phase factor in Ea. l26|) . being a topological quan- 
tity is closely related to the concept of the geometric 
Berry phases. To explicate this w we write explicitly the 
composite structure of the physical electron field using 
Eq.(|T3: 



Ca£(rr) = exp 



I I dr'VDiirr') 





emB(rr) j^^^^^y (28) 



The first term in the exponential in Ea. (|28|l is the usual 
dynamical phase factor where (j)e{rT) — yDf(rr) and 
(f)f,{v0) — (t>i{vO) . The second one, in turn, is the non- 
integrahle Berry phase factor:^^ 7£B(rT) = 2'KTmi{v) / (3 , 
where mp{r) marks the integer U(l) topological number. 



C. Topological ground state degeneracy 

It is known that the ground state degeneracy may arise 
from broken symmetries. Here, we argue that the ground 
state degeneracy of the charge states in strongly cor- 
related system states is a reflection of the topological 
properties of the system. The existence of topological 
features implies that the quantum eigenstates are not 
single- values under the continuation of the parameters in 
the Hamiltonian. Consider the "free" part of the action 
in Ea. (|21|l describing the dynamics of the U(l)gauge field 



dT- 



-E 



d(f>i{rT) 



dT 



2/i ^ 1 dcjjeirr) 
dT 



(29) 
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FIG. 2; Energy levels of the quantum Hamiltonian, Ea. (l3Ul . 
Half-odd integer values of the "statistical" parameter 
yield ground state degeneracy. 



This is the action of a particle moving in a plane around 
"magnetic flux " $/$o = 2/i/?7. The Hamiltonian corre- 
sponding to this action is simply 



d 



with the eigenenergies given by 
U 



E 



2^ 



2n 

T7 



(30) 



(31) 



The energy in Ea. (|31|l can be interpreted as the square 
of the kinetic angular momentum of a set of quantum ro- 
tors divided by the "moment of inertia" I — 2/U and the 
allowed angular momenta are uniformly displaced from 
integers by , which may be any real number. The 
energy spectrum is labeled by integers m£(r) and for inte- 
ger 2/i/[/ the ground state is non-degenerate. However, 
for the half-odd integer the ground state is dou- 

bly degenerate, see Fig|21 In this case the destructive 
interference between even and odd topological sectors is 
responsible for the ground state degeneracy. As we shall 
see in the following, this degeneracy will be crucial in ex- 
plaining the occurrence of superconductivity and anoma- 
lous properties on the verge of Mott transition, where the 
correspondence between the filling factor and the ground 
state degeneracy will also be established. 



quasiparticles which are characterized by effective elec- 
tronic parameters. Consequently, the tight-binding in- 
terpolation of the electronic structure is often used for 
fitting the experimental ARPES data for cuprates. In 
this type of analysis a conceptually simpler band theory 
is used to reveal how the strongly correlated electronic 
effects can be taken into account via the influence of the 
electronic band parameters. In the context of the present 
work in order to establish a link between the "bare" band 
parameters of the high-energy model in Eq.(^) and the 
"dressed" one of the low energy models we have to per- 
form the averaging over the fluctuation of the U(l) gauge 
phase fields according to 



[VjVf] 



-SleIIJ] 



,(32) 



[Vcfy] [VfVf] 
where the low energy action Sle is given by 

SleU, f] = -\n f [Vcb] e-5[^'/"-^]. (33) 



Performing the cumulant expansion in Ea. H33|) we can 
deduce in the lowest order the corresponding low energy 
fermionic Hamiltonian that has to be compared with the 
original one in Eq.QJ: 

+ E i'Vl,(r)/a.(r')- JE^1(^'^')^^(^'^') 



{{r,r')),a {rr') 

+ J2 t\{vv')flXv)U+,{v') 

(r,r') ,Q 



(34) 



where the dressed parameters encapsulating the effect of 
Coulomb interaction are given by 

tl(rr') =t^(rr')(e-''['^^(''^)~'^^+i('''^)l) (35) 
and (...) refers to the averaging over the U(l) phase field. 



V. BARE AND DRESSED BAND 
PARAMETERS 

With the help of the angle-resolved photoemission 
spectroscopy one gets a direct access to the density of 
low-energy electronic excited states in the momentum- 
energy space of Cu-0 planes in cuprates. Obviously, 
all the interactions of the electrons are encapsulated in 
ARPES data, but these are still difficult to evaluate. 
Since the underlying band structure of the hare electrons 
is a priori unknown, one way to think about these inter- 
actions is to consider simply the electronic excitations as 



with, the strongly fluctuating kinetic part (on the en- 
ergy scale U) So[<j)\ given by Eq.(|25J). On average, the 
effect of this renormalization due to the presence of the 
phase-phase correlation functions is the effective mass 
enhancement of the carriers as a result of the band nar- 
rowing, so that the "dressed" band parameters t*^ (where 
tx ~ t,t' ,t±) are used to match electronic spectra using 
the Zow-energy scale t—J model.^-^^ Typically, in cuprates 
t* ~ 0.5eV, t"'/t* ~ 0.15-0.35 and t*j_ is of order of mag- 
nitude smaller then the in-plane hopping parameters'^ 
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VI. PSEUDOGAP AND PHASE STIFFNESSES 
A. RVB pairs: single-particle "d-wave" gap 



where A(rr; r'r) is the complex pair field that is non- 
local in space. Furthermore, the corresponding effective 
action is of the form 



A routine Hubbard-Stratonovich decoupling is applied 
to the resonating valence bond (RVB) term in Eq. H25|l to 
give 



exp 



"'0 f /,.'\ 



e (rr') 

■VA*VA] I e-5[/-/.A*,A] 



(37) 



S[fJ,A\A] = f dr\A,{vr,v'r)\'-YY. f^ri 



A*(rT, r'r) 
V2 



[/T(rr)/Hr'r)~/x(rr)/T(r'T)] 



AeivT, r'r) 
V2 



[/T(rr)/^(r'r)-/x(rr)/T(rV)] 



(38) 



A saddle-point method can applied to the action in 
Eq.(j2Hl to give the self-consistency equation for the 
momentum-dependent gap parameter |A(k)| belonging 
to the Cu-0 plane. Assuming that A(k) is not chang- 
ing along the c-dircction, we can drop the layer index 
for this quantity. Out of the many possible mean-filed 
translationally invariant solutions in the RVB theory the 
"TT-fiux" phase is selected here because of its relation to 
the "d-wave" symmetry of the pseudogap in cuprates. It 
is governed by the equation 



-E 



ry2(k) 



N ^ 2E{k) 



tanh 



(39) 



where, 7y(k) = cos(fcj,a) — cos(fc,ya) with the quasiparticle 
spectrum of the fermionic part of the electron composite, 



and 



£;2(k) = [eUk)-/2p + |A(k)|^ 



A(k) = A[cos{akx) — cos{aky 



(40) 



(41) 



The gap parameter is a quantity with the short-range 
property, lim|r_r'|-,oo ^(r ~ r') = 0, essentially tied to 
local correlations on neighboring sites. As we see in 
the following, the presence of the "d-wave" pair function 
A(k) is not a signature of the superconducting state - it 
merely marks the region of non- vanishing phase stiffness. 
However, the presence of A(k) should be visible e.g., in 
ARPES spectra that picture the momentum-space occu- 
pation and therefore can detect the dispersion i?(k) with 
a gap for single particle excitationsiii 



B. Microscopic phases stiffnesses 

We can take advantage of the effective fermionic action 
in Eq. (|38() which is now quadratic in the Grassmann field 
variables that can be integrated out without any difficulty 
yielding a fermionic determinant: 



-So[0]+TrlnG-i 



(42) 



where 

G-^ =G-^ -T = {l-TGo^G-^ (43) 

is the effective propagator of the theory, while 
T^[T]„, (rr,rV) = 

+ t^{v - r')e-*[*^("^")-*^'('"'")l""M|^_r|4} az5{T - t') 
= f + f'+fi_ (44) 

is the matrix composed of the hopping integrals and 
phase factors. Furthermore, 



G- 



(rT,r'T') = 



d . \ , A,(rr,rV) . 



A*(rT, r'r) 

' 



5u,5{t~t') (45) 
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stands for the inverse of the "free" fermion propagator 
containing the gap A field. Here, 



2 ^ 
2 ^ 



(46) 



where CTq (a = 1, 2, 3) are the Pauh matrices acting in the 
Nambu spinor space, so that 



GoivTv'T') 



g(rrr'T') J^(rrr'T') 
^■^(rTr'T') ~g{r'T'rT) 



, (47) 



Using the self-consistency solution, Ea. l|39|l . and Fourier 
transforming to the frequency and momentum domain 
one obtains 



|A(k) 



Go(k.„) = ,(48) 

where Vn — {'^n + l)'"'//? are the (Fermi) Matsubara fre- 
quencies, n = 0,±1,±2, .... Now, expanding the trace 
of the logarithm in Eg. 142(1 we obtain up to the second 
order in the hopping matrix elements 

TrlnG-i = TrG'^ - ^Tr (GoT)^ - Tr (GoT') 



1. 



Tr GoT, 



(49) 



Finally, by performing summations over frequencies and 
momenta that are implicitly assumed in the trace oper- 
ation in Ea. l49|) we obtain an effective action expressed 
in the U(l) phase fields 



5, 



phi 



i0,^(rr) + |i<^,(rr) 

- ^ ^1 (A) cos [20, (rr) - 20,(r'T)] 

(rr'> 

- ^i(A)cos[0,(rr)-0,(r'r)] 

((rr'» 

- ^ :^_L ( A) cos [20, (rr) - 20,+i (rr)] | , (50) 

In our case U ^ J J a. so that the microscopic phase 
stiffnesses can be regarded as residual interactions to the 
dominant kinetic term of the phase model in Ea. l|50|) 
. This justifies the retention of only the lowest order 
non-vanishing terms in the electron hopping t and t' in 
Eq. (|49|l since for the ensuing microscopic phase stiffnesses 
one has, eg.-^^^^ ~ t/U and ^^-4^ ~ i' /U , respec- 
tively. All the stiffnesses in Eq.jOTj) (see also Ea. HA2p 
in Appendix A) rest on the single-particle gap due to 
the in-plane momentum space pairing among fermionic 
parts of the electron composites governed by the AF ex- 
change J: when A(k) = all the phase stiffnesses col- 
lapse. While J\\ and J±^ depend on the square of the 



corresponding hopping elements, the stiffness J^' is dif- 
ferent: it depends linearly on t' and governs the process of 
correlated particle-hole motion. Collective pair transfer 
events are costly for large J7, so that excitonic coherent 
charge transfer dominates the in-plane charge motion4^ 
The inter-plane stiffness is essential, however, in es- 
tablishing bulk superconductivity via the Josephson-like 
interplanar coupling. 



VII. FROM THE MOTT INSULATOR TO 
SUPERCONDUCTOR 

A. ODLRO vs. charge frustration 

Because of the composite nature of the electron field 
in strongly correlated system the occurrence of supercon- 
ductivity requires both the condensation of the of fermion 
pairs described by Ji[yt) , f]{vT) as well as the phase 
coherence which follows form the condensation of "flux 
tubes" e*'^*^'"'^-' attached to the /-fermions. Thus, non- 
vanishing of the pair-wave function A is not a sufficient 
signature of the off-diagonal long-range order (ODLRO). 
We can deduce this relationship from the definition of the 
superconducting order parameter which implies, 

^'«'(rT,r'r) ee (c4£(rr)c|£' (r'r')) 

= (5«,(/x,(rr)/T,(r'r')e-'[^^('-^)-^^('-'^')l) 
-^fe.A(rr,r'r)i^2 (5^) 

where V'o = (e"^*^'"'^)). The condensation of the "flux- 
tubes" from the electron composite has a transparent 
physical explanation. The phase factors which are intro- 
duced into the hopping elements by the gauge transfor- 
mation in Ea. H13|l frustrate the motion in the fermionic 
subsystem. However, when charge fluctuations become 
phase coherent^ as signalled by (e"^'^'"'^)) ^ 0, the frus- 
tration of the kinetic energy is released. The role of the 
gap parameter A also becomes apparent: pairing among 
fermionic parts of the electron composites is a necessary 
precondition for the existence of the microscopic phase 
stiffnesses (c.f. Ea. (|A2|) '). and, thereby, for the whole su- 
perconducting order. The opposite is obviously not true: 
the pseudogap state with A ^ may be phase incoherent. 
However, the appearance of bulk phase coherence in the 
presence of large Coulomb interaction U , whose energy 
scale by far exceeds that of microscopic phase stiffnesses, 
is not so obvious. Therefore in the following we elucidate 
the instrumental role of doping, since for the supercon- 
ducting order to occur the system should be brought in 
the vicinity of the degeneracy point, that is on the brink 
of change of the topological order. 



B. Effective NLcrM 

To proceed, we replace the phase degrees of freedom 
in Ea. H50|l by the unimodular complex scalar zi{rT) = 
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unity 



fields via suitable Fadeev-Popov resolution of 



1 ^ / [V'z]l[S{\ze{rr)\'-l) 

} 



(rr) 



(52) 



where the unimodularity constraint can be imposed using 
Dirac (5-functional, thus bringing the partition function 



into the following form: 



1 e 



(53) 



The (non-linear) unimodularity constraint can be conve- 
niently resolved with the help of a real Lagrange mul- 
tiplier A, so that the phase action in Ea. H5U|l can be 
suitably expressed by the effective non-linear cr-model 
(NLcrM) represented by the action 



S[z,z*] ^ r^M-E rdr'z?(rr)7-i(r-r')z£(rT)-2^|(A)^ 

£ -^0 \ ri (rr'> 



-zI{vt)zi{t't) -f C.C. 



((rr')> 



-zI{yt)zi{v't) -I- C.C. 



2:^±(A)^ 



-zl(rT)zt+i(rT) + C.C. 



(54) 



where 



variables becomes 



7(r - t') = — / [I?0]e'[^^("--)-^^('--')le-'^p"'[ 
Zq J 



(55) 



is the phase-phase correlation function calculated with 
the action in Ea. l29|l (see Appendix B). Since a part of 
the action in Eq. 154|) is quartic in the unimodular z-fields 
we employ the mean-field like decoupling 



-z^{rT)zi>,{v'T)+c.c 

{z*^{vT)zi,{v'T) [zI{tt)zi>{v't) + C.c] 



(56) 



to perform the closed-form integration over z-variables 
in Ea. H53|l . Here, the average (...) should be determined 
with the resulting effective action. To justify Ea. (|56|l . 
we observe that, formally, Eq.jS^Jl follows from the de- 
coupling of quartic terms in Eq.jMJ with the help of 
suitable Hubbard-Stratonovich transformation and sub- 
sequent use of the saddle-point method with respect to 
the emerging auxiliary variables. Next, we introduce the 
Fourier transformed variables 



1 



(3NN^ 



■^^z(q,w„)e" 



■i{u 



T,r — k-r— fcsc£ 



,(57) 



where q = (k, kz) with q — (k^, ky) and kz labelinq "in- 
plane" momenta and the wave vectors associated with 
the third dimension along c-axis, respectively. With the 
help of Ea. (|57|l the resulting quadratic action in the z- 



where 



r-l(q,c^„) = A- J(q)+7-l(co„) (59) 

is the inverse propagator for the z-fields and 

J{f\) = J]|(A) [cos(fcj,a) -f cos(fcya)] 
+ 2J^| (A) cos(A;a;a) cos{kya) 
+ J±{A)cos{kzc) (60) 

is the dispersion associated with the microscopic phase 
stiffnesses, where 

j;{A)=J(^{A){z*irT)zeir + l,t,r)) 

MA) = ^^(A)(z;(rr)z,+i(rr)). (61) 

Furthermore, Joi^n) is the Fourier transform of the bare 
phase propagator in Ea. H55|l . see Appendix B. 

C. SC critical boundary 

At the critical boundary demarcating the long-range 
ordered phase-coherent true superconducting state the 
static and uniform "order parameter" susceptibility di- 
verges, so that the suitable condition that can be read 
off from the action in Eq. H58|l is 



r-i(o,o)| 







(62) 
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which fixes the Lagrange parameter A at the transition 
boundary and within the ordered state. The parameter 
Ac is given by the solution of the uni-modularity con- 
straint equation: 



1 



1 



f]NN_i 



(63) 



The emerging ground-state phase diagram is depicted in 
FigEl It exhibits a periodic arrangement of phase in- 
coherent Mott-insulating lobes with the superconducting 
state above and between them. Apparently, at the de- 
generacy point defined by 



U 



(64) 



the superconducting state is most robust. Clearly, the 
above picture resembles that of a system described by 
the well known Bose-Hubbard HamiltonianS as a generic 



Hamiltonian for strongly correlated bosons. It covers the 
physics originating from the competition between the re- 
pulsive and kinetic term of the Hamiltonian, whose mag- 
nitude are proportional, in the present setting, to the pa- 
rameters U and phase stiffnesses, respectively. Moreover, 
this similarity is not accidental, since this is just an ex- 
ample of statistical transmutation where bosons emerge 
as fermions with attached "flux tubes" as a result of the 
gauge transformation^ 

Within the phase coherent superconducting state order 
parameter is given by 



1-1^2 



(3NNj 



(65) 



With the aid of Ea. HBl|l and IjBSp . by performing the 
summation over Bose Matsubara frequencies, we obtain 



1 - 



coth 



—y 



/?[/ 



/ 2[^(0)-^(q)] 



coth 



2^1 



T 



2m\ 



+ h 



2^' 



u 



u 



(66) 



where h{x—l/2) ~ x~ [x], while [x] is the greatest integer 
less than or equal to x. The remaining summation over 
the momenta can be efficiently performed by resorting 
to the lattice density of states as explained in Appendix 
C. The parameter ■0 as a function of temperature and 
the chemical potential is depicted in Fig01 which shows 
substantial enhancement of this quantity near the degen- 
eracy point. Note that the three-dimensional anisotropic 
lattice structure is essential, since even a very small in- 
terplanar coupling renders the phase transition in the 3D 
universality class as observed in cuprates^Thus, the ab- 
sence of t± will suppress the bulk critical temperature to 
zero, because for isolated stack of two-dimensional layers 
the NLcrM strictly predicts Tc = 0, in agreement with 
the Mermin- Wagner theorem. 



VIII. TOPOLOGICAL CRITICALITY AT THE 
DEGENERACY POINT 

In the preceding paragraphs we have shown that a 
theory of strongly interacting electrons can be trans- 
formed to an equivalent description of weakly interacting 
fermions which are coupled to the "fluxes" of the strongly 
fluctuating U(l) gauge field. In regard to the non- 
perturbative effects, we realized the presence of an ad- 



ditional parameter, the topological angle d/2TT = 2/i/C/, 
which related to the chemical potential. We argued that 
the configuration space for the phase field (p consists of 
distinct topological sectors, each characterized by an in- 
teger, entering the weight factors in the functional inte- 
gral and counts the topological excitations of the system. 
On the other hand, the ground state degeneracy depends 
also on the topology of the configurational space and the 
transition we encountered at 0/2Tr = 1/2 corresponds to 
an abrupt change of the ground state that is not related 
to any visible symmetry breaking. However, the exis- 
tence of different ground states that are related to the 
topological properties of the interacting electronic system 
is a hallmark of the topological order^^^ The latter is not 
associated with the symmetry breaking pattern, so it can- 
not be characterized by conventional order parameters in 
the Landau sense. We argue that the ground state de- 
generacy can be parametrized by a topological order pa- 
rameter being the average of the topological charge, i.e., 
the elements of the homotopy group of the U(l) gauge 
group and this parameter has a direct physical signifi- 
cance: in the large-C/ limit the electron density (i.e. the 
filing factor) is just given the mean topological charge 
rather than the number of fermionic oscillators. Further- 
more, in analogy to the Landau theory where the diver- 
gence of the order parameter susceptibility signals the 



0.50 



FIG. 3: The ground state phase diagram resulting from the 
effective action in Eg. lSUl . Here, in the filhng-control transi- 
tion, the control parameter is the chemical potential, which 
is conjugate to the carrier density. The picture shows the 
arrangement of Mott-insulating (incompressible) lobes Ml.m 
with topological order characterized but the winding number 
m — and m — 1, respectively with the phase coherent su- 
perconducting ground state between and above them. For 
large Coulomb energy U the phase coherent state is only pos- 
sible in the vicinity of the degeneracy point = 1/2. The 
curves are plotted for different ratios of the inter- to intra- 
layer couplings as input parameters: J±/J\\ = 0.001,0.01 
and 0.1 (from the top to the bottom) and show the prolif- 
eration of the superconducting state as the stack of coupled 
two-dimensional planes system crosses from 2D to 3D behav- 



FIG. 4: (Color online) Superconducting order parameter ip 
that signals the global phase stiffness shown as a function of 
the temperature and chemical potential, for U = 4eV, J = 



O.lSeV , t* = 0.5eV, t'* /t* = 0.25 and t\ 



■ O.OleV. 



number of of electrons per site in the Cu-0 plane) 

ne = ^Yl {cai{rT)cae{^T)) (67) 



consists of the fermion occupation coming from the 
fermionic part of the composite and a topological con- 
tribution resulting from the "flux-tube" attachment: 



phase transition between states with different symmetry, 
to indicate the change between different topologically or- 
dered states, we introduce the topological susceptibility 
being a derivative of the mean topological charge with re- 
spect to the statistical angle 9. Its divergence is related 
to the existence of distinct "vacua" , which cross in energy 
at the degeneracy point. Subsequently, we show that the 
topological susceptibility has a direct physical relevance, 
since it is related to the charge compressibility. It di- 
verges at the degeneracy point at T = and thus defines 
a novel type of topological quantum criticality, beyond 
the Landau paradigm of the symmetry breaking. 



A. Topological charge and the electron density 

In addition to the Coulomb energy U and temperature, 
the chemical potential /i plays a crucial role in Mott tran- 
sition, since it controls the electron filling Ue- An imme- 
diate implication of the composite nature of the electrons 
is that the electron occupation number (i.e., the average 



[^Cae{rT)cae{rT)j = /„^.(rT)/ai.(rT) 
2 /dM^T) 



iU \ di 



(68) 



The appearance of the topological contribution in 
Ea. (|68() is not surprising given the fact that "statistical 
angle", see Ea. H26|l . depends on the chemical potential 
and the occupation number is just its conjugate quan- 
tity. Owing that the U(l) topological charge (the wind- 
ing number) is given by 



TO^(r) 



1 



2-K 



27r 



dT4)t (rr) 







<pto(r)+2iTmt(r) 



d4>i{vT) (69) 



the mean value of the density of the topological charge 
can be written after performing Legendre transformation 
as 



2^l 2 /\d(l)t{YT) 



U U \ i OT 



(70) 
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Therefore the average electron occupation number Ue is 
given by 



Uf + nb — 



2/x 

77' 



(71) 



In the hmit of strong (weak) correlations rie interpolates 
between topological ni, (fermionic Uf) occupation num- 
bers. Clearly, in the large-C/ limit /i — > nfU/2, so that 
He — > fih and the system behaves as governed entirely by 
density of topological charge. The latter behaves in the 
large U limit as the typical density of hard-core bosons 
showing characteristic "staircase" behavior, see Figs. IS} 
13 Indeed, in this limit the system is described by the 
quantum rotor action in Ea. H29|l . in which the probability 
distribution function of the density of topological charge 
is gaussian and the problem has single-site character that 
can be analytically solved in a closed form: 
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(72) 



by making use of the Jacobi theta-function identity 
dJiiiv.q) _ ^ (-l)™47ri<7'^ 



i{iv,q) 



■sinh(27rmw). (73) 



The calculation of the mean topological density for the 
interacting problem, i.e. with the full phase action given 
by Eq. H50|l involving phase stiffnesses is a bit demanding 
since spatial correlations have to be included, as well. 
However, we can resort to the unimodular-filed NLcrM 
description given in Ea. (|52() . The result for rn, both 
within the Mott lobe and in the superconducting region 
is given by 



rit,(A), within MI 



within SC 



(74) 



where ip is the order parameter given by Ea.()66|l. while 



nh(A) = nobil-j) - 



_i_^|coth ^(^ 



/ 2mO)-^(q)] ^,^^^,/|X ,/| 



coth 



(75) 



with 5\ = A — Aq. Here, the parameter A is self- 
consistently determined via Ea. (|63|) whereas Aq is given 
by the solution of Eq. (|62|l . The summation over the wave 
vectors can be conveniently performed with the help of 
the lattice density of states. 



B. Topological susceptibility and charge 
compressibility 



Mott insulators have a clear distinction from metals by 
vanishing of the charge compressibility at zero tempera- 
ture, while this quantity have a finite value in metals^ In 
physical terms, the charge compressibility measures the 
stiffness to the twist of the phase of the wave function 
in the "imaginary time" direction. As we have shown, 
in the limit of strong correlations the physical proper- 
ties of the system are governed by the fluctuations of the 
topological charge. Thus, the effects connected with the 
nontrivial topological configurations of the gauge fields 
can be tested by performing the second derivative of the 
free energy with respect to the statistical parameter 0, 
see Eg . (1271) . which gives the topological susceptibility)^ 
i.e., the connected part of the two-point correlator of the 



topological charge densities at zero momentum: 

1 >. - d'^p{m) 
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(76) 



In Ea. H76|l . as a result of the non-trivial topology in 
the group manifold group caused by the non-simply con- 
nected structure, the partition functions Z{m) are given 
by the functional integrals taken over the field configu- 
rations in the topological class m only. The full parti- 
tion function Z, as defined by Ea. (|15|l involves all topo- 
logical sectors. Since the statistical angle parameter 9 
(and thereby the chemical potential /j,) acts as a ground 
state selector, the topological susceptibility can be con- 
veniently employed to detect transition between different 
topologically ordered states. Since these are labeled by 
the average topological charge, the abrupt change of this 
quantity will be signalled by the divergence of Xt at the 
degeneracy point 
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FIG. 5: Evolution of the fermionic occupation number n/ 
with increasing correlations at T = and for t* — leV, t* — 
0.5eV, t* = 0.25eV, and t* = 0.125eV from the upper left 
to the lower right. For all plots U = 4eV, J — 0.15eV and 
t'* /t* — 0.3. For large values of the Coulomb to band energy 
ratio the fermionic filling factor behaves as n/ ~ 



where /.tc is the value of the chemical potential at 
the degeneracy point. The topological susceptibility in 
Ea. (|76|l can be directly linked with the physical quanti- 
ties, namely the charge compressibility 



(78) 



which expresses the total density response of the sys- 
tem to a local change of the chemical potential. It is 
related to the shift of the electron chemical potential as 
a function of electron density which can be measured e.g. 
through the shifts of spectral features in photoemission 
spectra*^ Taking the derivative of Eq. (|71|l with respect 
to the chemical potential we obtain 
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U duf, 
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While the fermionic contribution dnj/d^ is regular, for 
the bosonic part dni, / dfi one gets 



(80) 



U drib 

Therefore in the large t/-limit the charge compressibility 
is entirely governed by the topological susceptibility and 
serves to distinguish that xt is zero in a Mott insulating'^^ 
region while it remains finite superfluid region and di- 
verges at the degeneracy point. 



FIG. 6: Evolution of the average topological number nj, with 
decreasing correlations at T = 0: U — 4eV, U = leV, U = 
0.5eV, and U = 0.2eV from the upper left to the lower right. 
For aU plots J = 0.15eV , t* = 0.5eV, t'* /t* = 0.25 and 
t*j_ — O.OleV. For large values of the Coulomb to band energy 
ratio the electronic filling factor behaves as ~ nj. The 
nearly linear dependence of rib as a function of the chemical 
potential near the degeneracy point — 0.5 signals the 

presence of the global phase stiffness, see Eg. 1741 . 



C. Electron mass enhancement at the degeneracy 
point 

Another remarkable aspect of the transition from one 
topologically ordered state to another is the great en- 
hancement of the effective mass m* of the electrons due 
to the collapse of electron kinetic energies due to the 
formation of the degenerate state at 2fi/U = 0.5. To 
estimate the change of m* we calculate 



where 



d\i\{k)/dklW=o 



Jl— ^g-''[0«(rT)-0£(r'r)]^ 



1 



\r-r\-- 



(81) 



(82) 



is the band renormalization factor, see Eq.JSHJ, where d 
stands for the lattice vector connecting nearest-neighbors 
sites on a two-dimensional lattice. Figure|Slillustrates the 
evolution of the effective mass as a function of tempera- 
ture in the vicinity of the degeneracy point. Interestingly, 
at 2fi/U = 0.5 which marks the "topological quantum 
critical point" the electronic matter in its charge aspect 




FIG. 7: (Color online) Average topological number ni, as a 
function of temperature and chemical potential forC/ = 4eV, 
J = O.lSeV , t* = 0.5eV, t'* /t* = 0.3 and tl = O.OleV. 

is very "soft" (see, FigQ making it very susceptible to 
transformation into alternative stable electronic configu- 
rations, namely to superconductivity, which we are going 
to analyze. 



IX. PHASE DIAGRAM FOR CUPRATES 

There has been a considerable amount of contro- 
versy regarding the observed pseudogap phenomena in 
cuprates.^^ One general class of theories views the pseu- 
dogap phase as resulting from performed pairs^S*^ The 
cuprates, however, are not in the strict Bose condensation 
("local pair") limit, since the photoemission still reveals 
the presence of a large Fermi surface. Furthermore, in 
the Bose limit, the chemical potential would actually be 
located beneath the bottom of the energy band, which is 
also not the case. The other class of scenarios (coined as 
"competing order" ) consider the pseudogap as not intrin- 
sically related to superconductivity, but rather proclaim 
it as competitive with superconductivity. Most of these 
proposals involve either a charge density wave^ or spin 
density waveS^, usually without long range order. How- 
ever, if there is a phase transition underlying pseudogap 
formation, a direct thermodynamic evidence (i.e. non- 
analytic behavior of the specific heat, the susceptibility, 
or some other correlation function of the system) must 
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FIG. 8: The ratio of the efTective-to-bare electron mass pa- 
rameter m*/me as a function of the chemical potential/i in 
the vicinity of the degeneracy point for U = 4eV, J = 0.15eV 
, t* = 0.5eV, t'*/t* = 0.25, t*j_ = O.OleV, and different tem- 
peratures T = IK, 50K, 115K and 300K from the top to the 
bottom. 



show up in existing experiments. Unfortunately, none 
of the spectroscopic data support a picture where the 
pseudogap phase represents a phase with true long range 
order. For the temperature-doping phase diagrams the 
two delineated above scenarios generally predict that the 
pseudogap characteristic temperature T*{x) merges with 
Tc{x) on the overdoped sideSS ("precursor" scenario) or 
T*(x) falls from a high value at low doping, compara- 
ble to the exchange energy to zero at a critical doping 
point inside the superconducting dome^ ( "competing or- 
der" picture). Below we show that both scenarios are 
consistently accommodated within the presence of topo- 
logical order, degeneracy point and accompanying phase 
coherence around it, as shown in Fig^| We sec the evo- 
lution of the charge compressibility k as a function of the 
chemical potential from the Mott insulator— with k = 
(at 2n/U = 1) to a point of degeneracy on the brink 
of the particle occupation change at 2/i/[/ =1/2 where 
K = cxD at T = 0. This is also the point on the phase 
diagram from which the superconducting lobe emanates. 
It is clear that, the nature of the divergence of k here has 
nothing to do with singular fluctuations due to sponta- 
neous symmetry breaking as in the "conventional" phase 
transition. Rather, this divergent response appears as a 
kind of topological protection built in the system against 
the small changes of /.t. Further, k — > oo implies that 
d^/due becomes vanishingly small at T = which re- 
sults in the chemical potential pinning, as observed in 
high-Te cupratesiSM^ 



A. Low energy scale pseudogap temperature T* 

In the pseudogap state at high temperatures one thus 
finds the coexistence of two distinct components; a 
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FIG. 9: (Color online) Contrasting behaviors at the degener- 
acy point T = 0, 2fi/U = ^ and around it. Right picture: the 
density plot of the electron effective mass parameter ml /me 
as a function of temperature and the chemical potential show- 
ing the strong depression of electron kinetic energies. Left pic- 
ture: The superconducting lobe with the density plot of the 
charge compressibility k diverging at the degeneracy point. 
In the "V-shaped" region fanning out to finite temperatures 
the electronic matter in its charge aspect is very "soft" ( i.e. 
highly compressible as opposed to the Mott state). The val- 
ues of the microscopic parameters for creation of the plots are 
the same as in Fig|3] 



state with gaped fermionic excitations (described by the 
fermionic part of the composite electron) and incoherent 
charge excitations (given by the attached "flux-tube"), 
which, as the temperature is lowered, enter the super- 
conducting state. The pseudogapped state is largely un- 
affected by the superconducting transition and does not 
participate directly in the superconducting behavior. As 
explained above the underlying mechanism for the ap- 
pearance of a gapped state with non- vanishing A is inti- 
mately connected with the antiferromagnetic correlations 



represented by the AF exchange J. We identify the tem- 
perature for which A sets in with the pseudogap tem- 
perature T* . This is also the temperature at which the 
microscopic phase stiffnesses in Ea. (IA2|) vanish, so that 
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Using results of Sec. VIA, the fermionic filling factor nj 
defined by Ea. H23l) can be computed explicitly as 
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With the help of Eq. (|83|) and using Eq. H39I) and we obtain 
that at T* 
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By eliminating the chemical potential from Ea. H85|) we 
get 
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revealing that T* is a universal function of the fermionic 
filing number Uf. Approaching half-filing [uf = 1) we 
can infer from Ea. (|86|l 



lim kBT*{nf) 
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41n(3) 



0.228 J; 



ksT^^axi (87) 



where T^^^ is the maximum value of the pseudogap tem- 
perature T*. For J = O.lSeV we obtain T^^^ = 396is:, 
see Figini 



Effect of doping on AF exchange and x 
phase diagram for cuprates 



T 



It is well known that the antiferromagnetic exchange 
J originates from the interplay between on-site repulsion 
([/) and the delocalization energy [t). The effect can be 
derived straightforwardly by expanding the energy to the 
second order in the hopping matrix element. This in- 
volves virtual double occupation and can be represented 
by an exchange process taking place on neighboring lat- 
tice sites: 




which yields a contribution ~ /U . There are also pro- 
cesses that are prohibited by the Pauli exclusion principle 
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such as I t,T) 0. However, with increasing hole dop- 
ing (i.e., departing from the half-fiUing) a given electron 
has fewer neighboring electrons to pair with, which re- 
sults in degradation of the exchange process described 
in Ea. (|55l) . Therefore, the AF exchange is leading to 
an effective interaction)S& which is a steadily decreasing 
function of x, vanishing at the critical doping Xc ■ 

J,«{x) = .7(1- Kx), (89) 

where K is treated as the lattice connectivity and the 
factor of K = 4 is the coordination number on the 2D 
square lattice. Therefore, the vanishing of JcS deter- 
mines the critical doping value Xc = 0.25. It is clear 
by inspecting Ea. H39(l that diminishing of the AF ex- 
change with doping will spoil the RVB pairing of the 
fermionic part of the electron composite by suppressing 
the c?-wave gap function A(k). Since phase stiffnesses in 
Ea. (|A2|l rest on A(k) the doping dependence of Jcff(2:) 
can be directly translated into the calculation of the x — T 
phase diagram that involves the doping effect on J^g , see 
FigEl By comparing Figllll and Fig^jwe can clearly 
see that the diminishing of the superconductivity in the 
overdoped region is just the result of the pair-breaking 
effect triggered by the doping dependent AF exchange. 
To summarize: T*{x) demarcates the region of non-zero 
microscopic phase stiffnesses which persist in the region 
characterized by the non- vanishing of the spin gap, as ob- 
served in high-frequency conductivity measurements.-^^ 
The origin of the spin gap is purely electronic and re- 
sults from the restricted space of available states that 
strongly correlated excitations on neighboring sites en- 
counter. It is described by the resonating valence bond 
singlet d-wave spin pairing of the fermionic part of the 
electron composite and is controlled by the AF super- 
exchange parameter J. The coincidence of the Fermi 
surface with the minimum gap locus as obtained from 
ARPES measurements^ also supports of a pairing gap 
interpretation of the pseudogapi^ 

C. Crossover to the "strange metal" state at Tj,: 
high energy scale feature 

The doping-dependent characteristic temperature T* 
in Ea. H8t)|) at which this pseudogap opens is in the un- 
derdoped region significantly larger than Tc. The phys- 
ical reason for this is transparent: T* marks the region 
of non- vanishing phase stiffness, albeit without global 
phase coherence (that appears at much lower temper- 
ature Tc). However, when the copper oxide supercon- 
ductor is driven in the normal state by applying a high 
magnetic field, a clear pseudogap feature at a similar en- 
ergy scale to the superconducting gap is observed in the 
quasiparticle tunnelling spectra and the pseudogap fea- 
ture persists up to the highest applied fields and does 
not depend on the magnetic field. Surprisingly, the Hall 
coefficient does not vary monotonically with doping but 
rather exhibits a sharp change at the optimal doping level 



for superconductivity™ This observation would support 
the idea that two competing ground states underlie the 
high-temperature superconducting phase. From this per- 
spective one has to conclude that any prospective order 
consistent with these observations implies that the pseu- 
dogap coexists with superconductivity and is essentially 
unchanged by a large applied external field. It is clear 
that sudden onset of the pseudogap at critical doping 
right at the point where the rigidity of the condensate 
wave function is at its maximum would be very diffi- 
cult to reconcile with the precursor scenario, but is very 
consistent with the onset of correlations which compete 
with superconductivity. This is precisely the outcome of 
the topological order which differentiates the electronic 
ground state into two states labeled by the topological 
winding number and with the degeneracy point sepa- 
rating them. It controls a remarkable concurrence be- 
tween normal state properties and the ground-state su- 
perconducting order setting up a unique critical doping 
point in the phase diagram where the transport prop- 
erties change very suddenly and where superconductiv- 
ity is most robust. Therefore we identify the crossover 
line Tg{x) where the charge compressibility undergoes a 
sudden change as an additional boundary hidden in the 
cuprate phase diagram, see Fig^l It is important to 
realize that in contrast to T*{x) the crossover line Tg{x) 
is controlled by the highest energy scale in the problem, 
namely the Mott scale set up by U . This explains why 
the anomalous behavior still persists even the supercon- 
ducting order is suppressed, eg. by the strong magnetic 
field. Considering the Hall effect one has to conclude 
that the mobile carrier density varies considerably with 
both doping and temperature. On the other hand, the 
ARPES Fermi surface remain essentially unchanged in 
the whole of the metallic doping range suggesting a con- 
stant density of states. These apparently contradictory 
results could be reconciled by observing that ARPES is 
sensitive to the momentum-space occupation and there- 
fore detect the excitation described by the fermionic part 
of the electron composite, whereas the charge transport 
properties are governed mainly by the "fiux tubes" which 
constitute charge collective variables. Given fact that the 
inverse of the Hall coefficient is proportional to the num- 
ber of carriers 1/Rh ~ "-e and Ue is governed by the 
topological charge ni, it is apparent why 1/Rh jumps in 
the vicinity of the QCP that is of topological origin. 



X. QUANTUM PROTECTORATE 

It is often difficult to formulate a fully consistent and 
adequate microscopic theory of complex cooperative phe- 
nomena and great advances in the solid-state physics 
are to a great extent due to the use of simplified and 
schematic model representations for the theoretical inter- 
pretation. In particular, the method of model Hamiltoni- 
ans has proved to be very effective. However, as it was re- 
cently argued^^, ab-initio computations have failed com- 
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FIG. 10: (Color online) Finite temperature phase diagram: 
the superconducting lobe Tc{(,) translated from the chemical 
potential = /i, upper right panel) to the particle occupation 
number = nj, lower left panel). Shaded area: the density 
plot of the charge compressibility k — Uk/2. The degeneracy 
point (T = 0, ^ = ^), where k, diverges transforms into the 
critical line on the nt — T phase diagram. Upper left panel: k 
and rii, as a function of temperature for T — O.IU showing the 
transition from the incompressible Mott state at 2fi/U = 1 
to a highly compressible region around the degeneracy point. 
The values of the parameters for creation of the plots are the 
same as in Fig|5| 



pletely to explain phenomenology of high Tc cuprates: 
it would appear that not only the deduction from micro- 
scopies failed to explain the wealth of crossover behaviors, 
but as a matter of principle it probably cannot explain 
it. Therefore it was concluded that a more appropriate 
starting point would be to focus on the results of ex- 
periments in the hope of identifying the corresponding 
"quantum protectorates" - stable states of matter whose 
generic low energy properties, insensitive to microscop- 
ies, are determined by a "higher organizing principle and 
nothing else" . From this perspective, each protectorate 
can be characterized by a small number of parameters, 
which can be determined experimentally, but which are, 
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FIG. 11: (Color online) Temperature-doping {x — 1 ~ rib) 
phase diagram calculated for the doping independent AF ex- 
change J (for the values of the parameters, see Fig|5Jl. De- 
picted is the superconducting lobe (bounded by the critical 
temperature Tc{x)) and the pseudogap temperature T* which 
marks the region of non-vanishing microscopic phase stiff- 
nesses (below T*(x)). Shaded area: intensity plot of the 
charge compressibility k diverging along the critical T = Q 
line (see, previous figure). 



in general, impossible to calculate from first principles]^ 
However, as we saw, a system with many microscopic 
degrees of freedom can have ground states whose degen- 
eracy is determined by the topology of the system. Pro- 
totypes of this kind of systems are provided by fractional 
quantum Hall effect. For example, the ground state de- 
generacy in FQH liquids is not a consequence of sym- 
metry of the Hamiltonian. It is robust against arbitrary 
perturbations, even impurities that break the symmetries 
in the Hamiltonian. Thus the topological ground state 
degeneracy on non-trivial manifolds provides a precise 
theoretical distinction between a topological and conven- 
tional order. The Hilbert space of quantum states de- 
composes into distinct topological sectors, each sector re- 
maining isolated under the action of local perturbations. 
This is a signature of its topological nature. Choosing 
the states from ground states in different sectors protects 
these states from unwanted mixing through the change 
of system parameters- protection within the sector is se- 
cured through a gapped excitation spectrum. In partic- 
ular, we found that for strong correlations the system is 
governed by the topological Chern numbers. However, 
the Chern number is a topologically conserved quantity 
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FIG. 12: Temperature-doping phase diagram as in Fig llll but 
with doping dependent AF exchange parameter Jeff (a;) given 
by Eg.dMl with K = 0.25. The temperature Tg{x) marks the 
sudden onset of greatly enhanced charge compressibility k in 
the "strange metal" state as departing from the half-filling 
{x = 0). Note, that the characteristic energy scale of these 
charge fluctuations (controlled by the Mott energy U ) is much 
larger than that of spin gap, controlled by AF exchange J. 



and is "protected" against the small changes of system 
parameters. Being an integer it can not change at aU if 
it has to change continuously. However, changing the in- 
teraction by a large amount may cause abrupt changes in 
ground state properties described by a different topolog- 
ical quantum number, which leads to a change of topo- 
logical order. This kind of stability might be generic for 
quantum systems governed by topologically non-trivial 
groups manifolds. Therefore one is left not only with 
the low-energy principle (the classic prototype being the 
Landau Fermi liquid), but the emergent physical phe- 
nomena are regulated also by topological principles that 
have a property of their insensitivity to microscopies and 
this quantum protectorate functions under certain topo- 
logical environments, through conserving of topological 
charges. 



XI. SUMMARY AND DISCUSSION 

In the present work focusing on the t — t' — t± — U— J 
model it is shown that the topological excitations of 
charge given by the collective U(l) phase field in a form 
of "flux tubes" attached to fermions can reproduce many 



robust features present in the phase diagram of high- 
Tc cuprates, thus substantiating one of the emerging 
paradigms in the condensed matter physics, namely the 
ubiquitous competitions in strongly correlated systems. 
The fundamental entities that carry charge (and spin) in 
the copper oxides are no longer the usual Landau quasi- 
particles but the "flux tube" fermion composites, so that 
the charge is no longer tied to the Fermi statistics. When 
charges are "liberated" then they can condense leading 
to superconductivity. This picture naturally leads to the 
pseudogap physics that is observed in the underdoped 
cuprates, which originates from the momentum pairing 
(in a d— wave pattern) of the fermionic parts of the elec- 
tron composite controlled by the antiferromagnetic su- 
perexchange J. This underlines the necessity of the fun- 
damental concept of fermion pairing in achieving the su- 
perconductivity. In the mathematical structure of the 
theory the gauge field is governed by the U(l) Chern- 
Simons term in the action of purely topological nature. 
From the canonical point of view the Hilbert space of a 
quantum theory has a non-trivial structure marked by 
the topological sectors which corresponds to a set of de- 
generate ground states. The topological ground state 
degeneracy provides a precise theoretical distinction be- 
tween a topological and conventional order: states with 
the same conventional broken symmetries may still be 
distinguished from each other on the basis of whether 
they are characterized by different topological quantum 
numbers captured by the homotopy theory gauge group 
manifold. In this paper, it has been shown that physics 
of the Mott transition is successfully covered in a topo- 
logical framework. The natural order parameter for the 
Mott transition is the topological charge related to elec- 
tron concentration for the filling-control scenario that 
selects topologically ordered states. Expectation value 
of the density of topological charge, determines also a 
dominant contribution to the topological susceptibility 
which, in turn, is related to the charge compressibility 
of the system. It diverges at the degeneracy point at 
zero-temperature and defines a novel type of topological 
quantum criticality, beyond the Landau paradigm of the 
symmetry breaking. Although the charge compressibility 
is completely suppressed in the Mott insulator because 
of the Mott gap, the criticality on the "strange metallic" 
side can be described adequately by the divergence of the 
charge compressibility at zero temperature. This critical 
enhancement of the density fluctuations extends to finite- 
temperature and is controlled by the Mott energy scale. 
This gives rise to another crossover line hidden in the 
cuprate phase diagram, where the charge compressibility 
undergoes a sudden change to the "strange metal" state. 
The crossover is governed by the Coulomb energy U, so 
that the density fluctuations at the instability towards 
the superconductivity surpass the effects of the spin fluc- 
tuation mechanisms governed by the antiferromagnetic 
exchange extensively studied for the cuprate supercon- 
ductors. This clearly demonstrates the inseparability of 
the high-energy Mott scale from the low-energy physics 
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in the cuprate problem and redefines the role of the chem- 
ical potential from a quantity that simply demarcates the 
boundary between filled and empty states to a selector 
of topologically ordered electronic ground state. 



APPENDIX A: MICROSCOPIC PHASE 
STIFFNESSES 

The microscopic phase stiffnesses to the lowest order 
in the hopping amplitudes are given by 
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Explicitly, after performing frequency and momentum 
sums in Ea. HAl|l we obtain 
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Here, e(x) — \J [i^ ^ jApx^ and 



p{x) = (l/7r2)K(v/l-(:rV4)) 



p{x) = - (2/7r2)E(Vl - Wl^)) (A3) 

where K(j:) and E(x) are the complete elliptic integrals 
of the first and second kind, respectivelyi^^ 



APPENDIX B: PHASE-PHASE CORRELATION 
FUNCTION 



where 1^3(2, q) is the Jacobi theta function, ^'^ which comes 
from the topological part of the functional integral over 
the phase variables. The function 193(2,(7) is defined by 
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and is /3-periodic in the imaginary-time r as well as in 
the variable 2/i/[/ with the period of unity. Fourier trans- 
forming one obtains 



m— — ( 

where 

Zq = exp 



1 +2^ # exp 



!3U 



4 



l-4[(m-f) 



(B3) 



(B4) 



is the partition function for the "free" rotor Hamiltonian 
in Eg. pijl . 

APPENDIX C: LATTICE DENSITY OF STATES 

In this Appendix we give the explicit formulas for 
the densities of states (DOS) for the anizotropic three- 
dimensional lattice that is helpful for evaluation of the 
sums over the momenta that appear in Section VII and 
VIII of the present paper. Our starting point is the 
dispersion relevant for the two-dimensional lattice with 
next-nearest interactions: 

E{V}) — — 2t cos(a/ca;) — 2t cos(&fcy) 

-f 4t' cos(afc;r) cos(aA:y) = 

— 2t [— cos{akx) — cos{aky) 

+ r cos{akx) cos{aky)] . (CI) 

The choice of such a dispersion is obviously motivated by 
its relevance as a simple means of modelling the quasi- 
particle band in the high-Tc cuprates. The density of 
states reads: 



piE) = 



TT / a 



TT / a 



-TT/a (27r/a) (27r/a) 



(C2) 



with 



By performing the functional integration over the 
phase variables in Ea. (|55|l we obtain 



7-i(r-r') = 



i9 



e W 



U 

xexp<; 



(Bl) 



P(£) 



K 



4-(e-r)3 
4(l+re) 



[9 (2 + r - e) X 



e (e + r) + e (-r - e) e (e + 2 - r)] , (C3) 



where e = E/2t and r — 2t' /t. Here <d{x) is the Heavy- 
side (unit-step) function. The expression in Ea. (|C3|l is 
valid for r < 1. For r > 1 one has to make the replace- 
ment p{e) |3fie/5(e)|. 
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The effect of inter-planar (c-axis) interaction can be 
incorporated as well via the following dispersion 



Performing the integration over we obtain 



EsdO^) = E{-k) - 2t, cos{cq,). 



(C4) 



In the presence of the system is a tree-dimensional 
anisotropic one for which the density of states becomes 



""'^ dk^ r'" dky r^" dq. 



-TT/a (27r/a) J-TT/b (27r/a) J-tt/c (Wc) 



(C5) 



/5(e)e[r2-(6-0^ 



(C6) 



in a form of the convolution involving previously calcu- 
lated DOS in En. deal. 
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